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I. INTRODUCTION 



The subject of integrable model is very fascinating largely because of its innumerable symmetries and a special class 
of solutions known as soliton solutions. Only a few systems in the nature are known to be integrable. The coupled 
| nonlinear Schrodinger equation (CNLS) and its higher order generalization namely coupled higher order nonlinear 
Schrodinger equation (CHNLS) are some of the examples of integrable equations which have direct relevance in 
the propagation of optical solitons in Kerr type nonlinear fibre. Coupled integrable systems have many important 
applications in photorefractive crystals jl| and also in all optical computations H. The cross phase modulation (CPM) 
phenomena in CNLS equation provides an interesting pulse shepherding effect to align the time arrival of the pulses 
||. The CPM phenomena along with group velocity dispersion (GVD) can also be utilized in compressing optical 
Dulses at the one soliton level 0. Recently the existence of multi-component solitons are experimentally established 



■ In this paper we have considered some interesting features of CHNLS equation, which describes the propagation of 
J> ' the optical pulses of very short width, of the order of lCP 15 sec. We have studied the novel inelastic collision of two 
t-H , solitons for an n-component CHNLS equation. The inelastic collisions among the three solitons are also studied. The 
three soliton interactions may be interpreted as a combination of two-soliton interactions occuring at three different 
points. This is in compliance with the three particle interactions in integrable systems. The energy exchange among 
the components of three solitons may occur at any of the three two-soliton interaction points. In this context, it is 
^v^j , important to note that the phenomena of shape changing vis a vis inelastic collisions and their consequences have 
been observed for the CNLS equation by considering two soliton solutions of both two and three components and 
subsequently n component generalization of two soliton solutions and n-complexes have been studied |2j. But the 
soliton solution being nonperturbative, the shape changing phenomena of CNLS equation does not ensure the same 
for CHNLS equation. Moreover, from integrability point of view study of the collisions of three solitons becomes 
indispensable and one needs to address this issue seperately. 



The paper is organized in the following sequence. In section I the CHNLS equation is introduced and its N 



soliton solution is given explicitly. In section III two-soliton interactions and its asymptotic analysis are obtained 
showing the novel inelastic collisions among the components of the solitons. Some interesting features of three soliton 
interactions, related to the integrability of the n-component CHNLS equation are addressed in section [V. Section |v| 
is the concluding one. 
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II. CHNLS EQUATION AND N-SOLITON SOLUTION 

The n-component CHNLS equation incorporating the effects of Kerr type nonlinearity and stimulated Raman 
scattering may be written as 

n n n 

E lz + %E lTT + 2i{J2 E*Ej)Ei + eE lTTT + 3e(^ E*Ej)E iT + 3e(^ E*E ]T )E, = (1) 

3=1 3=1 3=1 

where is n-component electric field, z and r denote the direction of propagation and time variable respectively. 
The parameter e is the ratio of the width of the spectra Auj to the carrier frequency u> such that e = — < 1 . 
A particular gauge equivalent form of (|l|) may be written as 0] 

n n 

qix + equtt + 3e^(q*qj)q it + 3e y^(g*gjt)gj = (2) 

3 = 1 3 = 1 

where the transformation relations are given by 

E t (z,T) = q l (x 1 ty(-^ x -^ 
X = z 

1 

t = T Z 

3e 

The transformed equation (||) is known as the coupled complex modified KdV equation (CCMKdV) whose N soliton 
solutions for n-component field is known . Notice that the CCMKdV equation in (0) is convenient particularly for 
the Lax representation of the system and consequently for applying the inverse scattering method to obtain ./V-soliton 
solution. However, as a result of the transformations (||), there is no change of the envelope function, but for a 
constant shift in velocity for all solitons. Consequently, the shape of the pulses for CCMKdV equation and CHNLS 
equation remains the same i.e., \Ei\ = |<&| and the transformation will not affect the subsequent results of the two 
gauge equivalent systems. 

The AT- soliton solution of the n component CCMKdV equation may be written in the most compact form as f| , 

N 

qi (x,i) = -2^(BC- 1 ) u e-^* t (3) 

3 = 1 

where, B and C are respectively n x N and N x N matrices whose explicit forms are given by 

(B) ii =iCW H (0)e- 8 ^ 3 -^ t (4) 

and 

(C) -VV CffilpWQ) e -»(V + A*-2A fc)W A3-A f)3; 

1 hj ^ h <+m + i( A 3*R'+in + i(A*) ' (A, - A*)(A, - A*) *° ( j 

In the above equation otij are the elements of the (n + 1) x (n + 1) scattering matrix and A is the spectral parameter. 
C\n+i)p 1S re l ate d to the elements of the scattering matrix at the position of the simple poles, wj^, = at(n+i)p(^|)U=o 
and / over a denotes the derivative with respect to A. 

The one soliton solution (1SS) for n-coupled system directly follows from (||,|],[5]) by considering N = 1 and is given 

as 

A 1 ) p-S^Lpimi 

Qi{x,t) = — — w — (6) 

cosh^ifl. + 



where rjm and rju respectively denote the real and imaginary parts of 771 = — 2i\\t — 8ieA* 3 x + and e fil1 = 
The complex variables arc defined as 



2 " A?~ 



n r (l) 
V^l U («+l)p 



= 1 ' u n+ln 



+i(Ai) 
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and 

An = Ai — A^ 

For convenience let us define Ai = Al -"+' AlJ ; where the subscripts R and I denote the real and imaginary parts and 
consequently r)i R and r\\i become 

mi = -Aifli + eXiRCiXjj - X\ R )x + - (7) 

Vib = -Ai/[t - e(X 2 u - 3X 2 1R )x] (8) 

From the argument of the cosh function of the 1SS (||) and (||), it is straight forward to identify the width of the 
soliton pulse as T — | Axj | 1 and the soliton travels with a group velocity V g — [e(Afj — 3X\ R )]~ 1 in the positive x 
direction when X 2 j > 3X 2 R and in the negative x direction when X 2 j < 3X 2 R . 



III. TWO-SOLITON INTERACTION 



The two soliton solution (2SS) can be obtained from by putting N = 2, the explicit form of the i lh component 

of 2SS being 



2J2 2 -i (C ( ™ ) 1 -e , ' m + e rn+ ' n2+r >™ +& ™^j 



Det[C] 

where, 



(9) 



r,s— 1 



In ( |l0| ) Xij = Xi — X*_^ and r\ m = —2iX* m t — 8ieX^x + i^. Det[C] is the determinant of 2 x 2 dimensional matrix 
obtained from (C)y rtq). The form of Det[C] in terms of the spectral parameters emerges rather lengthy but it will be 
useful for the asymptotic analysis of the interacting solitons and will be seen later. In terms of the spectral parameters 
Det[C] is given as 



Det[C] = 1 + e Vr+vt+Rr S + e 2 mR +2 mn +Q ( H ) 

with 



r,s— 1 



e R - = (12) 

A rs 

n r (s) r *(r) 
_ V"^ ^n+lp^n+lp ^ 

p=l a '(r l +l)(n+l)( A *) a '(n+l)(n+l)( A -5) 

and = Det[/t]Z)et[£], [k] = (j 12 -) (no sum over [£] — (Ay) -1 , [k] and [£] being 2x2 matrices. To understand 
the interaction in a more explicit manner we analyse the asymptotic limits of the 2SS (^ . The asymptotic limit may 
be obtained by observing the 2SS when both the solitons are infinitely apart. This may be achieved by taking the 
limit r]2R — > ±oo. As a consequence i th component of the remaining soliton acquires the form 

= A ^' ±)e "* J (14) 

where A^^ and $( z± ) may be interpreted as the amplitudes and phases of i th component of I th soliton respectively. 
Notice that in general Af ' are different from A,[ . This may occur due to exchange of energy between the solitons. 
The amplitude Af ' in terms of C„l u is given as 

CO 



= ^ (15) 
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The expression for A^ + ^ however, is more involved and may be written in terms of Ay as 



A {l+) = A (J-) T i (16) 



T\ may be interpreted as the transition matrix and is defined as 




and 



y* y* y^* y^ 



(1-^A 2 ) 



1? ~ Vl-AiA 2 \Z ( S^7 )( a|^ ) (18) 
where Ai = Kl2 '^ 11 , A2 — K21 \ 22 , If |T?| = 1, the solitons pass through each other without being affected in their 

^11A12 ' ^22^21 III L o O 

shapes and sizes. Otherwise, we will see that the solitons exchange energy at the time of interaction. The phase shift 
3>W as a result of collision may be obtained from the following relation. 

$(0 = _ $(*-) (19) 

where = |(-Rn) and <I> ( ' +) = \{Q - R 22 ) and consequently the phase shift becomes $ w = |(Q - i?n - #22)- 

It is now interesting to note that for equal values oiC^\ u for each component of the solitons, \T-\ becomes unity and 



consequently there is no exchange of energy among the components of each soliton. However, if the relative phases 
among c£| H are introduced the energy of each individual component no longer remains constant due to collisions 
and in that case \Tj\ ^ 1. This analysis has been presented graphically by plotting two soliton interaction for a two 
component system with the s = 0.1, Ai = —1.5 — i and Xi = — .5 — 2i. One of the solitons is moving with group 
velocity —1.74 and the other is moving with a group velocity 3.07. In /inures ([!]) two soliton interaction is plotted 

with C^ +li — 1 for each component, where the energy profile of each component of the solitons remains unchanged. 

In /igwres^) two soliton interaction is plotted with nontrivial relative phases, C^}\ 1 = C^+i 2 = ^n+i 2 = ^ ano - 

^n+i 1 = 46(1 — i) showing a considerable exchange of energy among the components. /i<7Mre(||) shows that the overall 
intensity profile, represented by \q\ of the 2SS remains unchanged although there are appreciable energy exchange 
among the components of a soliton. It is an interesting result in the context of all optical computations |2j leading to 
the construction of the logical binary gates. 
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FIG. 2: Plot of |gi| and \q%\ showing inelastic collision 
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FIG. 3: Plot of \q\ showing constant energy profile 



IV. THREE-SOLITON INTERACTIONS 



The explicit form of the three soliton solution (3SS) can be obtained from (||,||||) by putting N = 3, the explicit 
form of the i lh component of 3SS being 



Det[C] 



3 



n+lp c 



V" 2r] iR + r)„ 

3 3 f~Lx 

(adj^) mm (adjjC. j ) mm e 2r > iR +^ + £ Det[k m ]Det[C m ]e %™ 



i,i,m=li^j^„ 



where, 



3 3 



£>ei[C] =-1- e^ +r >h R »-J2 E (ad?'K)ifc(ad?'>C) jfc e 2 *«+^ + ^ - £ ( a d7K) fefe (adjX) fefe e 2 "-+ 2 "^ 



i=i fcj=i 



i,j,fc— i 



+DetiK]Det[£}e 2m n +2n2 « +2v3 n 
with £ is a 3 x 3 dimensional matrix with elements (A" 1 ) and [k] is also a 3 x 3 dimensional matrix with elements 

3 J 

(no sum over [%?] is the matrix [«], where the jth row has been replaced by S e ji^n+ip an d 'C' 7 is the matrix 

where the jth row has been replaced by a unit row vector, X) e^. The phenomena of shape changing due to exchange 
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FIG. 4: Contour plot of three soliton interaction 



of energy among the components of solitons becomes more interesting in the three-soliton interaction case. It is 
observed that the interactions take place at three different space time points with two solitons interact at each time. 
Thus the energy exchage among the components in three solitons interaction may occur at three different points. This 
is demonstrated in /igureM) by the contour plot of the three solitons interaction with e = 0.1, Ai = —1.5 — i, A2 = 

-0.01 + 2., A3 = -0.9 + 2i,c|2 M =C<3i i=4+i 1 = 10- 3 ,C.^i 2=4+1 2 = 1 andCgij, = 10 3 -0.8i Two of the 
solitons are moving in the positive x direction with group velocities 2.5 and 6.36 and the third soliton is moving in the 
negative x direction with group velocity —1.74. The nature of three soliton interactions opens up another possibility 
of three solitons interaction where second and third solitons interact first instead of first and second solitons leaving 
the final configuration same in both the cases. 

V. CONCLUSION 

In conclusion, we have shown the novel shape changing phenomena asscoiated with two soliton interaction for n 
component CHNLS equation. The three solitons interaction also exhibits energy exchange, but the exchange may 
occur at three points. This may lead to a more flexibility in constructing logical gates in all optical computing systems. 
The three solitons interaction also demonstartes an interesting feature conforming the exact integrability of the system 
a la Zamolodchikov and indicates the existence of Yang Baxter like relation, which will be publishedelsewhere. 
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